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We calculate radiative corrections to the masses of the Higgs bosons in a minimal supersymmetric 
model that contains an additional non-anomalous (7(1)' gauge symmetry. With some fine-tuning of 
the U(l) charges of the Higgs fields, it is possible to suppress the Z-Z' mixing. We use this fact, 
along with the lower bound on the lightest Higgs mass after LEPII era, as a criterion to restrict the 
set of parameters in our analysis. We calculate the mass of the lightest Higgs and its mixing with 
the other Higgs bosons, in a large region of the parameter space. 



PACS numbers: 



I. INTRODUCTION 



m 
o 
o. 

Supersymmetry (SUSY) is one of the most viable extensions of the Standard Model (SM) for curing the quantum 
instability of the Higgs sector. SUSY, however, has a new hierarchy problem concerning the natural scale of the 
Higgsino Dirac mass - the fi problem. Many supersymmetric extensions of SM have been proposed to cure the p, 
problem. Although there is no natural solution to this problem in the minimal supersymmetric models, the next- 
10 \ to-minimal SUSY Q solves the fj, problem, but at the expense of inducing large-tension domain walls, which can 
over-close the universe. Among all possible extensions of supersymmetric SM, however, the ones that invoke an extra 
U(l)' symmetry, which forbids a bare /i term are particularly promising since problems with the generation of 
tadpoles and axion are automatically avoided. Also, extra U(l)' gauge symmetries, with their associated extra Z' 
bosons, naturally arise in effective theories coming from breaking of grand unified theories (GUTs) or in superstring 
, compactifications 3]. The above two reasons are the major motivations for our paper. 

In general, the presence of a Z' boson in the low-energy spectrum generates additional neutral current transitions, 
which can have important implications for precision tests Q or flavor violation [f| . It has been shown that the existing 
data allow for a light Z' boson with family-dependent couplings [||. In a general Z' model, the hierarchy of the soft 
masses in the Higgs sector leads to different vacua in which the Z 1 boson can be heavy or light, although the Z-Z 1 
mixing is sufficiently small in all cases. The analysis of Erler and Langacker Q shows that the Z' boson can be as 
^ 1 light as 200 GeV with a Z-Z 1 mixing angle of the O(10" 3 ). 
Qh, Although a low-scale extra U(l)' symmetry stabilizes the \i parameter to TeV scale, the tree-level Higgs potential 
is far from representing the physical observables with sufficient precision. In fact, experience from the minimal model 
shows that there are large corrections to the Higgs masses, once the radiative effects are taken into account J?]. This is 
• also the case with the Z' models, as the partial analysis of the Higgs sector shows 8]. Therefore, for a proper analysis 
of the collider data concerning the Higgs production, it is necessary to compute the Higgs masses and couplings 
•rH , at least at the one-loop level. In fact, this is necessary even for making comparisons or putting constraints on the 
/\ • parameter space using the LEP results @- 

Our main goal in this paper is to work out the neutral Higgs boson masses and mixing angles at the one-loop level 
in supergravity models with an additional U(l)' gauge symmetry. A general tree-level analysis of such models can 
be found in .101. We will assume that the parameters of the model have already been fixed to the TeV scale via 
the one-loop renormalization group equation (RGE) running from the string scale with appropriate initial conditions. 
The low energy mass spectrum and appropriate vacua that naturally suppress the Z-Z' mixing angle have already 
been derived in ^(J f° r the tree-level Higgs potential. We will further assume that there is no violation of the CP 
invariance in the Higgs sector, which is a simplifying assumption rather than a result following from high-scale model 
building. 

In the next section we review the properties of the tree-level potential. In Sec. 3, we compute the effective potential 
and derive the Higgs mass spectrum at the one-loop level. We will also briefly comment on the relevance of two-loop 
effects. In Sec. 4, we present some numerical estimates of the masses and mixing angles for likely values of the 
parameter values. Here, the lower bound set on the lightest Higgs mass by LEPII is used as a constraint to restrict 
the values of our parameters. Of course, to improve the bounds on our parameter space and to further restrict the 
parameter values, we need to consider the Higgs production and decay rates. The dominant channels are the Higgs- 
strahlung (Bjorken) process, e~e + — > Z* — > Zjhk, and the Z decay into two Higgses, e~e + — > Z* — > hjhk- However, 
we will work out the details of these processes and cross-sections in a separate work. For now, the bounds that we 
obtain using the lightest Higgs mass as a constraint are still valid. Finally, in Sec. 5, we conclude the work and discuss 
its implications for Higgs phenomenology. 
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II. TREE-LEVEL EFFECTIVE POTENTIAL 



We first describe the structure of the effective potential at the tree level. The gauge group is the same as that of 
the Standard Model, but with an additional U(l)' factor, i.e., G = SU(3) C x SU(2) L x U(l) Y x U(1) Y ', with coupling 
constants g 3 , g 2 , <7i, and g[, respectively. The Higgs sector contains two Higgs doublets Hi and H 2 , and one singlet 
S. The matter multiplets are given by left-handed chiral superfields. For specific charge assignments of the chiral 
superfields with respect to the gauge group, we refer the readers to Ref. Hi) . For our purposes, we only need to know 
the charges of the Higgs fields under the extra U(l)' gauge symmetry. We denote the Higgs charges by Qi, Q 2 , and 
Q s . The superpotential is given by 



W = h s S(H 1 -H 2 ) + h t U c (Q-H 2 ) 



(1) 



The most important feature of the above superpotential that distinguishes it from either MSSM or NMSSM, is the 
absence of a cubic term in S and a term proportional to Hi ■ H 2 , usually called the fi term 16]. The gauge invariance 
of the superpotential under the U(l)' forbids the appearance of such terms. Although a /i term is absent from the 
superpotential, an effective /x parameter is generated by the vacuum expectation value of the scalar field S. We use 
the hatted fields to denote chiral Higgs superfields, and the un-hatted fields to denote scalar Higgs fields. 
We parameterize the explicit soft breaking of supcrsymmetry by 
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where Xi are the gaugino fields, and the hermitian conjugate terms are assumed to keep the potential real, 
tree-level Higgs potential follows from Csb, F and D terms: 
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where Qi are the U(l)' charges of Hi and Q s + Qi + Q 2 — due to the gauge invariance of the superpotential under 
the extra gauge symmetry. Above we notice that g[ always appears in combination with Qi. Hence, it is convenient 
to absorb it in the definition of Qi and define new charges Qi = g[Qi- Therefore, g[ will not explicitly appear in our 
formulae unless stated otherwise. 

We assume that all the coupling constants in the above potential are real. At the tree level, the potential cannot 
violate CP symmetry cither explicitly or spontaneously. A possible phase could come from Ah s ; however, such a phase 
could be absorbed into the global phases of the Higgs fields. At the one-loop level, CP symmetry can be explicitly 
broken due to the complex phases in the scalar quark sector. However, we set all the CP-violating phases equal to 
zero and consider only the CP-conserving scenario. 

The Higgs sector of the theory contains ten real degrees of freedom. Each Higgs doublet contains four real fields 
and the Higgs singlet contains two real fields. After the electroweak symmetry breaking, four of the ten fields become 
the longitudinal components of the four vector bosons in our model. The remaining six fields result in three scalars, 
one pseudoscalar, and one charged Higgs. We decompose the Higgs fields as 
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where the neutral components will be further separated into scalar and pseudoscalar bosons below. 
The vacuum state of the theory is defined by the Higgs vacuum expectation values (VEVs): \cj>\) 



(4>+) = 0, 



(^4>i) = Vi/y/2, and (S) — V3/y2, where Vi are real, and v 2 = vf + v 2 — (246 GeV) 2 . The effective \i parameter is 
generated by the VEV of S, and is defined by /i s = h s v^/^/2. Here Vi = wcos/3, v 2 = vsm/3, and v 3 = x. For this 
to be a physical minimum, the potential must be negative when evaluated at the point [vi ,v 2 ,x), and the masses of 
the Higgs bosons must be positive. Even when these conditions are satisfied, the above point is not guaranteed to be 
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the absolute minimum. Whether it is still acceptable depends on the location and depth of the other minima and the 
width between them. At the minimum point, the potential has vanishing first derivatives with respect to the three 
CP-even scalars, i.e., all tadpoles vanish. This enables one to trade the soft mass-squared parameters for their VEVs. 

The tree-level masses of the Higgs bosons are obtained by diagonalizing their corresponding field-dependent mass- 
squared matrices. To this end, we need to substitute <j>® = (vt + <j)%+ itpi) /y/2 and S = (x + 4>3 + 2^3) /v2 into the 
potential. Here 4>% and tp% stand for CP-even and CP-odd directions, respectively. Using the basic definitions 
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we form the Higgs mass-squared matrix. Evaluation of Eq. Q at the tree level is straight forward. Defining m 2 = A/i s , 
we obtain 
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For the pseudoscalar mass-squared matrix, we get 
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The eigenvalues of the scalar mass-squared matrix correspond to the masses of the Higgs bosons. Although these 
eigenvalues can be obtained analytically, they are often too complicated to be useful. However, from the structure of 
the above matrix we can obtain useful information about its smallest eigenvalue, which corresponds to the lightest 
Higgs mass. Namely, for any symmetric n x n matrix, its smallest eigenvalue is less than the smaller eigenvalue of it 
left upper 2x2 sub-matrix. With this observation, we get 
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where Qh — Qi cos 2 (3 + Q 2 sin 2 (3. The first two terms are familiar from NMSSM f], while the third term is unique to 
the model under consideration. Notice that this term allows the lightest Higgs mass to be larger than that predicted 
by either MSSM or NMSSM. 

After appropriate rotations, the pseudoscalar mass-squared matrix gives one non-zero eigenvalue corresponding to 
the physical pseudoscalar mass. The other two eigenvalues which are zero, are the Goldstone degrees of freedom. After 
electroweak symmetry breaking, these become the longitudinal components of Z and Z' . The Z-Z' mass-squared 
matrix is given by 
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where M 2 , = (Q 2 v\ + Q 2 v 2 + Q 2 x 2 ) and A 2 = (1/2)^5? 
together with the Z-Z' mixing angle, are given by 



9l 



Qi — v\Q 2 \. The eigenvalues of the above matrix, 
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The mixing angle otz-z 1 has to be smaller than a few times 10 3 , so that Mz 1 would correspond to the observed Z 
boson mass. For completeness, we also give the expressions for the pseudoscalar mass and the charged Higgs mass: 
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sin p cos p 2 

From the above, it is clear that the pseudoscalar mass is never negative, while the charged Higgs mass can be lower 
than the W boson mass, and can even run to negative values for some choices of the parameters. 

In the next section, we include the main one-loop contributions to the tree-level effective potential. In general, 
the tree-level potential is written in terms of the running coupling constants and masses, which are defined at some 
renormalization point Q. However, the tree-level effective potential written in terms of running parameters is too 
sensitive to the choice of Q, and one cannot make reliable calculations. The situation is considerably improved when 
one includes the one- loop contributions to the effective potential We take into account the one-loop top/stop 
and sbottom effects, which are the main corrections. 

III. ONE-LOOP EFFECTIVE POTENTIAL 

As explained at the end of the previous section, the most important one-loop contribution to the tree-level effective 
potential comes from the top and scalar top quarks. However, the contribution of the bottom scalar quarks can also 
be sizable when tan/3 ~ 40 or larger. We take both contributions into account. For the rest of this section, we will 
state our results in full generality, making no assumptions about the numerical values of our parameters. 

The stop and sbottom mass-squared matrices are given by 
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where k\ = (g 2 - g 2 /3) /4 + Q Q Q U fc* = h\ - (g\ - g 2 /3) /4 + QqQ 2 , k\ = QqQ s , l\ = gj/3 + QuQi, it 
h 2 t - <?f/3 + QuQi, It = QuQs, and k\ = h 2 - (g 2 + g 2 /3) /4 + Q Q Q U k\ = (g 2 + gj/3) /4 + Q Q Q 2 , k b s = 



= ' l b — 5i/6 + QdQi, 1% = .9i/6 + QdQi, and finally, l b — QdQs- The cancellation of triangle anomalies gives 
Qq = ~ Qi/3, Qu = (Qi — 2Q2) /3, and Qd = (Qi + 2Q2) /3. The eigenvalues of the above matrix are the masses of 
the left-handed and right-handed stops and sbottoms, given by 

™L = \ trM i T \\j(trM q ) 2 -±detM q . (12) 

Using the stop and sbottom masses from above, we can express the one-loop correction to the effective potential by 
the Coleman- Weinberg formula [I3 

V 1 = k 




(13) 

where Q is the renormalization scale in the MS scheme and k = 3/(327r 2 ). We sum over q = (t,b) and j = (1,2). 
Here, rh 2 = h 2 \H$\ and to 2 = h 2 \H°\ . The one-loop scalar and pseudoscalar mass-squared matrices are given by 

(A4) =(^- §ij ^ 
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where the second term in the brackets is due to the fact that the position of the minimum has shifted because of the 
one-loop effects. By substituting Eq. (|12|) into Eq. I|13|) and substituting the resulting expression into Eq. (|14|> . we 

get 
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where q = t (b) for top (bottom) quark/squarks. Above, A 22 = 2hfv2i)2, X\i — 2hfviV\, and the rest of the A^- are 
zero. The values in parentheses correspond to the bottom quark/squarks. The k\ and Zf are defined above, following 
Eq. (|ll(l . The functions J 7 and <7 that appear in Eq. (|15(l . are the usual loop amplitudes which also appear in the 
MSSM effective potential. These are given by 
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where one particularly notices that T depends explicitly on the renormalization scale. By combining Eqs. JBJl, i|l(J|l . 
and 1|14[1 . we can diagonalize the total pseudoscalar mass-squared matrix to obtain the mass of the pseudoscalar: 
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(18) 



We can check the validity of our expressions for Mj> and M.ij by comparing them to the well-known MSSM results. 
Our model reduces to MSSM if we fix \i s — h s x/ v2 = /i and set h s = g[ = Q, = 0. In this limiting case, we identically 
recover the usual MSSM results computed in Q- We can also obtain a useful upper bound for the lightest Higgs mass 
at the one-loop level. 
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Although we did not take into account the two-loop effects, we do not expect these corrections to be very large. 
Indeed, as was shown in Ref. pj| for MSSM, the two-loop effects are less than a few GeVs. In our future work, we 
will include the two-loop effects in the model under consideration. But given the fact the not even the one-loop effects 
have been worked out in this model, it is important to have these corrections first before higher order effects are taken 
into account. 



We numerically diagonalize the total one- loop scalar mass-squared matrix, given by Eq. I|15fl. to obtain the mass 
of the lightest Higgs scalar. In order to get concrete results, we must fix some of the parameters in our model, which 
include (72, 9\, ht, h s , mjj, mo, itiq, At, Ab, A, v, x, tan/3, Q\ and Qi . The U{1)' gauge symmetry does not affect 
the mass of the W boson. This fixes g 2 = 0.65 and v = 246 GeV. According to the one- loop RGE analysis of [Icj . 
h t ~ 1.0 — 1.2 and 0.35 < h s < 0.9 — 1.0. We leave the soft supersymmetry breaking masses rriQ : u,D and A, A t ,b, h s , 
tan /3, and Q as free parameters against which the lightest Higgs mass will be plotted. Furthermore, because we are 
considering the electroweak symmetry breaking driven by a large VEV of the singlet field S, much larger than v, we 



The values of g' x , Q\ and Q2 are not directly constrained by the experimental data. However, in most GUT- 



motivated models with an extra U{\)' factor, g[ — ^^/5A 5 /3J g\ tan^, where g\ is the hypercharge gauge coupling, 

and X g is of the order of one with the exact value depending on how the GUT gauge group is broken down to the 
SM gauge group ^ij. For our calculations, we take 175. = 0-3- To fix the charges Qi, we take the Z — Z' mixing angle 
Q'z-Z' to be smaller than a few times 10 -3 . Without excessive fine-tuning of the theory, this requires h s ~ g[Qs 
(see [lOj for details). This together with 0.35 < h s < 0.9 imply Q s ~ O(l). Furthermore, the gauge invariance of 
the superpotential under the U(l)' gives Qi + Q 2 = —Q s - Because Q1Q2 > as is argued in [l(|, we conclude that 
\Qi\ ~ 1. A natural choice would be to take Q\ = Q 2 = — 1. 

However, according to Eq. 0, this specific choice restricts the value of tan/3 severely. In fact, az-z' < 3 x 10 -3 
implies 1 < tan/3 < 1.3. It turns out that tan/3 is highly sensitive to the ratio of the charges. Fortunately, the masses 
of the Higgs bosons are not very sensitive to the charges at all. This is because in the mass-squared matrix, each 
charge Qi is accompanied by a factor of g[. At the one-loop level, it is the top quark Yukawa coupling that gives the 
largest contribution to the Higgs masses. Therefore, we are free to choose Qi without effecting our results significantly. 
We have also verified this numerically. Because tan (3 is the only quantity that is really sensitive to the ratio of the 
charges, we choose charges that allow tan/3 to vary over a wide range of values. If we take Q\ — SQ 2 = —3/2, then 
we can have 1 < tan/3 < 60 without violating the bound on otz-z 1 ■ If we choose other values for the charges that are 
still of the order of one, the masses of the Higgs bosons would not change by a significant amount. 

We find that the Higgs masses are sensitive to the value of h s . For this reason, we calculate the Higgs masses for 
h s = 0.4, 0.6, and 0.8. Finally, we fix the scale at which we carry out the numerical calculations. We take the scale, 
denoted by Q not to be confused with the charges Qi, to be of the order of the electroweak scale. This is indeed 
necessary for the consistency of our analysis. We have also verified that as Q varies between 100 — 2000 GeV, the 
Higgs masses change by less than 2 — 3 GeV almost independently of the other parameters, see Fig. Etc). It is clear 
that the change in the lightest Higgs mass is because of different choices of h s , and not because of Q. That the Higgs 
masses are relatively stable against the choice of the scale, is a restatement of the stability of the effective potential 
against the scale when one- loop effects are taken into account, as we mentioned above. In the actual calculations we 
fix Q = 300 GeV. 

In summary, we are studying the lightest Higgs mass, which is phenomenologically the most interesting one, as 
we vary mu,D,Qi A ti b, A, tan/3, and h s . We have summarized our results in Figs. 1 and 2. We have plotted the 
lightest Higgs mass for 100 < A < 2000 GeV, -1600 < A t ,b < 1600 GeV, 200 < mu,D,Q < 1000 GeV, 1 < tan/3 < 60, 
and < h s < 1.0. We take mu,D,Q > 200 GeV to keep squark masses positive. A negative mass for even one of 
the squarks induces color breaking minima in the potential, which is clearly undesired. The above bounds provide 
a fairly comprehensive region for our parameters. The physically allowed values must respect the LEPII constraint, 
M hl > 114 GeV Q. 

Using Figs. ^and[21 we summarize our main results as follows: First, h s > 0.8 is excluded if tan/3 > 2, see Fig. 
|2Ie). For tan/3 < 2, h s can be 0.8 or larger. However, according to Figs. [IJa-f), when h s > 0.8, tan/3 > 4 is excluded 
unless A > 1300 — 1350 GeV. We conclude that h 8 > 0.8 and tan/3 > 4 cannot be allowed simultaneously. Second, 
the lightest Higgs mass, , is relatively stable against the soft supersymmetry breaking masses mjj,D,Q and against 
Ab, see Fig. EJa,b,c,e). There is some variation with respect to toq, to_d, and Ab, but only when tan/3 is very large. 
Mf ll varies quite a bit when A t and A vary. However, when tan /3 is not near its extreme values, either too large or 
close to 1, then Mh 1 is relatively stable against A t and A as well, see Fig. dd,f). Mh ± stays pretty much fixed as Q 



IV. NUMERICAL EXAMPLES 



fix x = 1.5 TeV. 




7 



varies between 200 — 2000 GeV. This is demonstrated in Fig. |2c), for different values of h s . We have also plotted 
the lightest Higgs mass as a function of fi s , Figs. |2Ja,b). Notice the extreme fine-tuning of the charges Qi and Q 2 
in order to have tan /3 large with no restriction on x, and respecting the experimental bound on the Z — Z' mixing 
angle. With Fig. |2Ib), we arrive at the same conclusion as above. Namely, for large tan/3, h s < 0.8. Larger values of 
tan /3, however, are not generally favored by the above model. 

We have verified numerically that for nominal values of tan/3 and h s , say, tan/? ~ 4 — 15 and h s ~ 0.4 — 0.6, the 
lightest Higgs mass is about 130 GeV, which is very close to the MSSM prediction. This can be seen from our graphs 
by examining the dashed black and red curves. Based on the MSSM analysis [13J, we do not expect the two- loop 
corrections to add more than a few GeVs to this value. That there is an agreement between the above model and 
MSSM is no surprise, because the above model is just an extension of MSSM. This is in fact a check on the validity 
of our results. 

The above conclusions might be slightly effected when CP-violating effects are taken into account. However, the 
inclusion of CP- violating phases coming from h s A s and h t A t mixes the scalar and pseudoscalar mass-squared matrices. 
The resulting mass-squared matrix is a 4 x 4 matrix giving four Higgs bosons with no definite CP properties. This 
requires a different analysis than what we have done in this paper. We consider CP violation in a separate work. 

Finally, for phenomenological reasons, it is desirable to know the Higgs mixing angles. In Fig. 2(f), we have plotted 
the three Higgs mixing angles against tan/3. The derivation of the Higgs mixing angles is shown in appendix A, 
which is already well known. From the graph, one can see that for smaller values of tan /3, the lightest Higgs state is 
dominated by a mixture of H® and , while its mixing with the singlet stays less than about 10 percent. However, 
when tan/3 is larger than about 10, the lightest state is dominated by H® , and the total mixing with and S stays 
less than 10 percent. This result is not surprising because larger tan/3 implies a larger VEV for H%, 

V. CONCLUSIONS 

In this work, we studied the one-loop effects on the lightest Higgs mass in a minimal supersymmetric model 
augmented by an Abelian U{1)' gauge symmetry. We calculated the top and stop/sbottom one-loop effects in the 
framework of the effective potential approach. The most important issue concerning the minimal models extended 
by a U(l)' factor, is the mixing of the Standard Model Z boson with the Z' boson associated with the extra gauge 
symmetry. In order for such models to be phcnomenologically viable, the Z—Z' mixing angle has to be very small, 
less than a few times 10 -3 . We used the smallness of the Z-Z' mixing angle, together with the lower bound set on 
the lightest Higgs mass, 114 GeV, by LEPII data @, to constrain our parameter space. We showed numerically that 
the one-loop effects due to the top and stop/sbottom quarks are non-negligible. 

The radiative corrections to the Higgs boson masses and mixing angles are crucial for interpreting and predicting 
the Higgs production and decay rates in upcoming colliders. In linear colliders, for example NLC, the main production 
mechanisms are the Bjorkcn process and pair-production process, each of which requires a precise knowledge of Higgs 
boson masses and their couplings to the gauge bosons. (For the analysis of these processes at the tree level, see e.g 
[l5|L The model at hand predicts a larger upper bound on the Higgs boson masses than MSSM. Therefore, even 
if the MSSM bounds are violated in the near-future colliders, the model at hand, which generates the [i parameter 
dynamically, will accommodate larger Higgs masses. 
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APPENDIX A: DERIVATION OF THE HIGGS MASSES AND MIXING ANGLES 

The Higgs mass-squared matrix is given by a symmetric 3x3 matrix, (Mij). The eigenvalues are given by the 
solutions of the following characteristic equation: a; 3 + rx 2 + sx + t — 0, where 

r = -tr(Mij), 

a = - l -[tr 2 {M l0 )-tr(Ml)}, 

t = -det(Mij). (Al) 
With the help of the auxiliary parameters p = s — r 2 /3, q = 2r 3 /27 — rs/3 + t, d = p 3 /27 + q 2 /4:, and i9 ~ 



arccos ( — qj y'—4p 3 /27) , we can express the Higgs masses as follows: 



1 j~p (d 2?r 



M, n = --r + 2W- £1 cos 



3 V 3 V3 3 
1 r~p ( 'd 2?r 



Mf l2 = --r + 2W-^cos 



3 V 3 V3 3 



where we require d, r, t < and s > 0, to ensure that the masses are physical, i.e., positive. Because the Higgs 
mass-squared matrix is real and symmetric, it can be diagonalized by means of an orthogonal transformation 0, 
where 

O t M vj O = diag (M 2 3 , M 2 2 , M 2 hi ) . (A3) 

In the basis (<j>i,4>2,(f>3) — {H^H^S), the mass eigensta mtes are defined by — Ylj^i^ij^ji f° r * = 1, 2, 3. 
More specifically, we have 

hi = 031^1 + 03202 + 033^3, 
h 2 = 02101 + 02202 + 02303, 

^3 = 01101 +01202 + 01303- (A4) 
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We are defining our fields such that < M% 2 < M% . Due to orthogonality, we have 0\ { + 0\ + 0^ = 1, for 
i = 1,2,3. Since we are mainly interested in the mass of the lightest Higgs boson, we only plot the values of 3 i, 
032, and 033, as these are the only mixing angles that determine the composition of the lightest Higgs boson. 
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FIG. 1: Black, red, and blue curves correspond to h s = 0.4, 0.6, and 0.8, respectively. 
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M, [GeV] vs. n [GeV] 
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FIG. 2: Black, red, and blue curves correspond to h s = 0.4, 0.6, and 0.8, respectively. 



